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ABSTRACT 

We use the Massachusetts Institute of Technology general circulation 
model (GCM) dynamical core, in conjunction with a Newtonian relaxation 
scheme that relaxes to a gray, analytical solution of the radiative transfer 
equation, to simulate a tidally locked, synchronously orbiting super-Earth 
exoplanet. This hypothetical exoplanet is simulated under the following 
main assumptions: (1) the size, mass, and orbital characteristics of GJ 1214b 
(Charbonneau et al., 2009), (2) a greenhouse-gas dominated atmosphere, (3), 
the gas properties of water vapor, and (4) a surface. We have performed a 
parameter sweep over global mean surface pressure (0.1, 1, 10, and 100 bar) 
and global mean surface albedo (0.1, 0.4, and 0.7). Given assumption (1) 
above, the period of rotation of this exoplanet is 1.58 Earth-days, which we 
classify as the rapidly rotating regime. Our parameter sweep differs from 
Heng and Vogt (2011), who performed their study in the slowly rotating 
regime and using Held and Suarez (1994) thermal forcing. This type of ther- 
mal forcing is a prescribed function, not related to any radiative transfer, 
used to benchmark Earth's atmosphere. An equatorial, westerly, superrotat- 
ing jet is a robust feature in our GCM results. This equatorial jet is westerly 
at all longitudes. At high latitudes, the flow is easterly. The zonal winds do 
show a change with global mean surface pressure. As global mean surface 
pressure increases, the speed of the equatorial jet decreases between 9 and 15 
hours local time (substellar point is located at 12 hours local time). The lat- 
itudinal extent of the equatorial jet increases on the nightside. For the two 
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greatest initial surface pressure cases, an increasingly westerly component 
of flow develops at middle to high latitudes between 11 and 18 hours local 
time. On the nightside, the easterly flow in the midlatitudes also increases 
in speed as global mean surface pressure increases. Furthermore, the zonal 
wind speed in the equatorial and midlatitude jets decreases with increasing 
surface albedo. Also, the latitudinal width of the equatorial jet decreases as 
surface albedo increases. 



Keywords: ATMOSPHERES, DYNAMICS ; EXTRASOLAR PLANETS ; 
ATMOSPHERES, STRUCTURE 
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1. Introduction 



Infrared observations of transiting planets over the past decade have 
open ed a window on a new regime of atmospheric dynamics (e.g. 



Knutson et al 



20071 ). The first transiting exoplanets discovered and modeled were "hot 
Jupiters", planets broadly similar in size to Jupiter that orbit their parent 
stars in close-in orbits (often with periods less th an a few Earth-day s ) anc 



hence receive an immense amount of stellar flux ( IMayor and Queloz 



1995 



Charbonneau et al 



2000l ). Their m asses and rad ii sugg est that their at 



mospheres are hydrogen-dominated (IFortney et al. 



20071 ). similar to Solar 



System gas giants. They are also expected to be tidally locked and syn- 
chronously orbiting their parent star, always keeping one face to their parent 
star. These properties imply that their circulation patterns will be entirely 
unlike anything found in the Solar System. 

1.1. Review of exoplanet general circulation models 

Many exoplanet general circulation mo dels (GCMs) have been develo ped, 



almost all from existing terrestrial GCMs. 



Showman and Guillot 



(120021) de 



veloped the first true 3-D GCM of an exoplanet (in this case hot Jupiters 
such as HD 20 9458b) using t h e Exp licit Planetary Isentropic Coordinate 



(EPIC) model (bowling et al 



19981 ). This model differed from previous 



models of brown dwarfs in that they included the intense irradiation re- 
ceived from their parent star and differed from previous exoplanet models 
in that they included the effects of advection. They predicted superrotat- 



5 



ing (westerly) winds, which later models have verifie d, and phase-curve off- 



sets b efore the latter was discovered observationally. ICooper and Showman 



(120051) modeled HP 2094 5 8b us ing version 2 of the ARIES/GEOS dynamical 



core (ISuarez and Takacsl . 



19951 ) with Newtonian rela xation to a prescribed 



equili 



Drium temperature. This model improved on IShowman and Guillot 



(120021 ) by having a more realistic radiative-equilibrium temperature, better 



resolution, and a deeper domai n. It was thus able to 



night temperature difference . 



Cooper and Showman 



j etter p redict the day- 



Cooper and Showman 



(120061 ) then used the 



(120051 ) model to study the chemistry of HD 209458b. 
CO and CH4 were represented as passive tracers. They found that where 
there is disequilibrium chemistry of CO and CH 4 , their abundance also de- 
pends on the meteorology of the planet. 

Several hot Jupiter models have used the equivalent barotropic formu- 
lation or shallow water version of the primitive equations of atmospheric 



fluids. 



ChoetaL 



( 120031 ) developed a GCM for HD 209458b using the equiv- 



alent barotropic formulation of the shallow-water equations with Newtonian 



heating to a prescribed 



Showman and Guillot 



equi 



ibrium state and obtained different results than 



((2002|). Namely, they predicted a ban ded structure 



with t hree broad zonal jets and easterly flow at the equat or. 



( 120031 ). in a companion paper, generalized the 



Cho et al. 



Menou et al. 



(I2003f ) GCM to 



other giant planets by varying the Rossby and Burger numbers. They found 
that close-in giant exoplanets lie in a different regime of flow patterns than So- 
lar System gas giant planets. Namely, the former have few zonal jets, while 
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the la tter have many. 



ChoetaL 



(120081 ) extended the study of 



ChoetaL 



(120031 ) to a much larger parameter space by varying the initial RMS velocity, 
global mean temperature at cloud tops, thermal forcing amplitude, Rossby 
deformation radius at the pole, and nondimensional Rhines length at the 
equator. They concluded that giant exoplanets have a polar vortex that 
revolves around each pole, a low (two or three) number of zonal jets (in con- 
trast with Solar System gas gian t s) , an d a temperature fiel d that depends on 



the circulation. 



Rauscher et al. 



(|2008|) used the GCM of 



to pre dict observational signatures of hot Jup iters 



Cho et al 



(2008) 



Langton and Laughlin 



Cho et al. 



(120031 ) and found 



(120071 ) used a shallow water model similar to 
a cold spot centered 76° east of the antistellar point. They attributed this 
disagreement to a difference in their radiative time constant. 

Meanwhile, hot Jupiter model s with other types of dy namical cores con- 



tinued to develop and improve. 



Showman et al. 



(120081 ) simulated the at- 



mospheres of the hot Jupiters HD 209458b and HD 189733b using a GCM 
based on the Massachusetts Institute of Technology (MIT) GCM dynami 



cal co re, which uses the full 3-D privative equations. The 



Showman et al 



(120081 ) GCM was radia tively forced using New tonian relaxation to a radia- 



tive equilibrium state. 



Showman et al. 



radiative transfer (RT) scheme (IMarley and McKay 



(200 9J) continued t his work using a 



19991 ) that calculated 



the RT fluxes directly using a correlated-k scheme with multiple wavelengths 
in the two-stream approximation. Both of these studies showed flow from the 
substellar to antistellar point along the equator and over the poles at lower 
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pressures, while a westerly equatorial jet and e asterly polar flow developed at 



higher pressures, 
dynamical core of 



Hoskins and Simmons 



Menou and Rauscher (2009) used the Intermediate GCM 



(119751 ). which solves the primitive 



equations. This GCM was thermally forced by Newtonian relaxation to a 
prescribed equilibrium temperature and used t o compare the barotrop i c (hot 



Jupiter) and baroclini c (Earth-like) regimes . 



used the same model as 



Rauscher and Menou 



Menou and Rauscherl (120091 ) t o check that hot Jupiters 



(120 loi ) 



were being mod eled consistently with the re s ults of 



2005 



2006) and 



Showman et al 



(12008 



20091 ) 



Cooper and Showman 



Perna et al. 



( I2010al ) used the 



Rauscher and Menou! (120101 ) GCM to show that magnetic drag is a pla usible 



me chanism to limit wind speed s in hot Jupiters 



the 



Rauscher and Meno 



u| ( 201oj 



Perna et al. 



feoiObl ) used 



GCM to show that Ohmic dissipati o n is a 



non-negligible heat source in hot Ju piters. The 



Menou and Rauscherl (120091 ) 



Rauscher and Menoul (120121 ) to include 



GCM was subsequently improved by 
dual band, double-gray RT. 

Several studies have investigated the effects of particular aspects of hot 



upite r GCMs on the resulting circulation and temperature. 



Burkert et al 



( 120051 ) used an inviscid, primitive equation model in 2-D (varying in azimuth 



and radius) with flux limited radiative diffusion to show that the nightside 



tempe rature is sensitive to atmospheric opacity. Similarly, 



Dobbs-Dixon and Lin 



( 120081 ) used a flux limited radiative inviscid hydrodynamical model (trun- 
cated at ±70° latitude) to show that atmospheres with large opacity have a 
large day-night temperature difference, while atmospheres with reduced opac- 
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ity have a more uniform temperature. 



Dobbs-Dixon et al. 



on the model of 



(120101 ) improved 



Dobbs-Dixon and Linl (120081 ) by using decoupled radiative 



and thermal components with updated opacities. They also added viscosity, 
and found that high viscosity resulted in subsonic flow, while low viscosity 



produced supersonic flow. The a ddition of viscosity changed t 



patterns and variability of flows. 



le circulation 



Langton and Laughlinl (120081 ) used a shal- 



low water model with two a two-wavelength radiative transfer to investigate 
the effect of high orbital eccentricity on hot Jupiters. They found that the 
atmospheric response was driven primarily by the intense irradiation at peri- 
astron and that the resulting expansion of heated air produced high velocity 



turbulent flow ( with some planets deye 



as well). Finally, 



oping superrotating acoustic fronts 



Thrast arson and Choi (120101 ). using a GCM with the prim- 



itive equations and Newtonian relaxation to a prescribed equilibrium state, 
found that different initial states of the initial zonal jet (which had the most 
significant effect), thermal drag timescale, spectral truncation number, and 
superviscosity co efficient, lead to mark e dly d ifferent results for circulation 
and temperature. iThrast arson and Choi (1201 ll ). using the same GCM, found 



that a short Newtonian relaxation time led to large amounts of unphysical, 
gird-scale oscillations that contaminated the flow field. 

Heng, Menou, and Phillipps (2011) performed benchmark simulations of 
Earth, tidally locked synchronously orbiting Earth, a shallow hot Jupiter, 
and a deep HD 209458b. The purpose of benchmark simulations is to com- 
pare dynamical cores of GCMs to determine if differences are due to nu- 
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merical modeling procedures or "higher level" physical parameterizations. 
Heng, Menou, and Phi llipps (2011) followed the standard Earth method of 



Held and Suarezl (119941 ) with a slight modification to include the effects of a 
tidally locked synchronously orbiting planet. This scheme used Newtonian 
relaxation to relax to the following radiative equilibrium temperature: 



T, 



eq,hs 



max < T< 



strat i 



T sur f — Aft cos (A — 180°) cos0 — A z log f — j cos 2 ^ 



p_ 



where T eq ^ s is the 



(1) 



Held and Suarezl (119941 ) equilibrium temperature (in K), 



T«trat = 200 K is the stratospheric temperature, T sur f=315 K is the surface 
temperature, = 60 K is the equator to pole temperature difference, A 
is longitude, (ft is latitude, A z = 10 K is a characteristic temperature dif- 
ference in altitude, p is pressure, p s is the surface pressure (here taken to 
be 10 5 Pa), R is the sp ecific gas constant, an d c p is the specific heat at 



constant pressure. Since 



Held and Suarezl (119941 ) were comparing terrestrial 



dynamical cores, they set their parameters in Eq. [T]to mimic the temperature 
distribution in Earth's atmosphere. Heng, Menou, and Phillipps (2011), in 
their benchmark test, compared the Geophysical Fluid Dynamics Laboratory 
(GFDL) - Princeton flexible modeling system (FMS) (see references in Heng, 
Menou, and Phillipps 2011) spectral and finite difference cores. They found 
qualitative and quantitative agreement between the two cores in all but the 
deep HD 209458b case, where horizontal dissipation changed the results with 
errors of 10% or greater in the temperature and wind fields. Heng, Frierson, 
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and Phillipps (2011) continued the development of the GFDL-Princeton FMS 
exoplanet model by adding dual-band radiative transfer, dry convective ad- 
justment, and a uniform slab of finite thickness to model the surface. Their 
goal was to bridge the purely dynamical 3-D simula tions with those th at in- 
corporate multiband radiative transfer. In addition, iPerna et al.l ( 120121 ) used 
the Heng, Menou, and Phillipps (2011) and Heng, Frierson, and Phillipps 
(2011) models to show that on hot J up iters, irradiation by the parent star is 
the main driver in radiative and dynamical properties, but that atmospheric 
opacity also plays an important role. 

Thus far in our discussion, we have focused on Jupiter-size planets or 
larger, but models for smaller planets exist too. The main commonality for 
all of these models is that the planet is tidally locked with and synchronously 
orbiting its parent star. There seems to be a continuum in size for these 
planets, at least in the modeling sense, ranging from Jupiter size or larger, to 
Neptune size, to mini-Neptunes, to super-Earths (two to ten Earth masses), 
to Earth or smaller. Moreover, there is a also a range of modeling choices 
on how to specify the bottom boundary condition (in terms of the heat and 
momentum exchange), which will be noted through out Section 



the 



urning our attention to smaller-sized planets, 



Lewis et al. 



Showman et al. 



( 12010h used 



(120091 ) GCM to study the hot Neptune GJ 436b and 
found that as atmospheric metallicity was increased, the circulation changed 
from one dominated by mid-latitude jets and minimal longitudinal temper- 



ature differences, to one dominated by an equatorial jet and large day-night 
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temperature differences. 



Menoul (12012ft used the lRauscher and Menoul (J2012J) 



GCM to model super-Earth GJ 1214b with H2O, solar, and x30 super solar 
metallicity and found that superrotating equatorial jets were robust features 
for these cases. 

Long before the discovery of super-Earth exoplanets 
used a terrestrial planet GCM (i.e. sh allow atmosphere cove ring; a thick rocky 



Joshi et al 



(1997) 



surfac e) , based on the simple GCM of 



Joshi et al. 



(|l995|) and 



James and Gray 



( 119861 ) to simulate a tidally locked, synchronously rotating exoplanet. They 
used a two stream, gray RT representing a gas with CO2 or H 2 radiative 
properties. The model was somewhat of a blend of terrestrial (radius, gravi- 
tational acceleration at the surface, gas constant, and stellar insolation) and 
Mars constants (for instance, the atmosphere could freeze and sublimate ac- 
cording to CO2 parameters), but also had a slow, Titan-like rotation rate 
(16 Earth-days). Their results contained a longitudinal, thermally direct 
cell, which transported heat from the dayside to the nightside. The mass 
was returned to t he dayside through the polar regions at low levels. In the 
Joshi et al.l (119971 1 model, the circulation aloft was fairly longitudinally sym- 



metric, and a superrotating jet was present at the equator. The winds were 



we sterly a 



most everywhere. 



Joshil (120031 ) followed up the work of 



the Intermediate GCM described in 



Joshi et ai- 



de F. Forster et al. 



1997) by switching to 



(120001 ) and adding 



more model features such as non-gray radiative transfer; clouds; continents 
and oceans; and evaporation, precipitation, and transport of water. They 
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found that the presence of a dry continent caused higher temperatures on 
the dayside and allowed accumulation of snow on the nightside. Furthermore, 
the absence of any ocean l ed to higher day-night temperature differences. 



Merlis and Schneider! (120101 ) used an idealized, terrestrial aq uaplanet 



primitive equation, spectral GCM (jO'Gorman and Schneider 



20081 ) that con- 



tained an active hydrologic cycle. They used a two stream RT scheme that 
had a shortwave component and two gray long wave components (one where 
optical depth depended linearly on pressure to represent a well mixed gas 
such as C0 2 , and the other with a quartic dependence representing H 2 0). 
They performed two simulations, assuming the planet was tidally locked to 
and synchronously orbiting its parent star: one in which the rotation period 
was equal to one current Earth year and the other in which the rotation pe- 
riod was equal to one current Earth day. They found that waves and eddies 
shaped the winds, temperature, and precipitation in the rapidly rotating at- 
mosphere, while simple divergent circulations dominated the slowly rotating 
atmosphere. 

Heng and Vogt (2011) used the Heng, Menou and Phillipps (2011) GCM 
to simulate super-Earth Gliese 581g. They varied the radiative cooling time, 
the rotational frequency of the planet, the mass of the planet, the density of 
the planet, the equator to pole temperature difference, the global mean sur- 
face pressure, and presence or lack of a troposphere. They found the radiative 
cooling time to have the most substantial effect on the results. A long radia- 
tive time scale compared with the advection time scale allowed the circulation 
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to diffuse the temperature field. 



references in 



Selsis et al. 



Selsis et all (120111 ) used the LMD GCM (see 



20111 ) to simulate a tidally locked, sychronously 
orbiting super-Earth sized planet with a CO2 atmosphere to determine the 
wavelength and phase changes in the thermal emission as atmospheric pres- 
sure increased. 

A fe w studies have i nvest igated the stability of liquid water on super- 



Earths. 



Pierrehumbert 



(1201 ll ) modeled Gliese 581g with the FOAM GCM 



(previously used to study sn owball Earth and oth er deep-time paleoclimate 



problems; see references in 



Pierrehumbert 



201lh . They found that their 



"Eyeball Earth" simulations, i.e. mostly frozen planets with a pool of liquid 
water centered at the substellar point c ould maintain liquid H 2 stability 



at the surface. 



Wordsworth et al. 



(120111 ) simulated a dif ferent super-Earth 



Gleise 581d, using the LMD Mars model (see references in 



Wordsworth et al. 



201 ll ) to determine if a CO2 atmosphere was stable against atmospheric col- 
lapse and if liquid water could be stable on the surface. They found this 
is the case for surface pressures with over 10 6 Pa of CO2 for both land and 
ocean surfaces and varying amounts of background gas. 

1.2. Present Study 

The motivation f or our work came wi t h the tr ansit observati o n of s uper- 



Earth GJ 1214b by 



Charbonneau et al. 



(j2009h 



Zalucha et al 



(120 lOh and 



Zalucha and Gulbisl (120121 ) had previously adapted the MIT atmosphere GCM 



to Pluto and Mars, respectively. Given the success of the Mars and Pluto 
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versions, which we call collectively the MIT small body GCM, it seemed rea- 
sonable to extend these GCMs to the super-Eart h regime with so l id or l iquid 
surfaces. Note that this is a separate effort from 
use the giant planet version that is derived f rom the MI 
R ecent spectral observat i ons of GJ 12 



Showman et al 



(120091 ) who 



2011 



Crossfield et al. 



2011 



Bert a et al. 



4b (|Bean et al. 



ocean GCM. 



2010 



2011 



Desert et al. 



20121 ) indicate th at the atmospheric 



composition is either metal rich or hydrogen-dominated (ICroll et al. 



2011 



de Mooij et al. 



20121 ) . H 2 and CO2 have been suggested as po ssible candi 



dates for the metal rich case. However, 



Rogers and Seagerl ( 120101 ) have shown 



using planet interior structure models, that if GJ 1214b is composed of wa- 
ter, it must include a super-critical fluid layer. Such a lay er is not currently 



presen t in the MIT small body GCM. A separate study by 



Nettelmann et al. 



( 1201 ll ) favored a Neptune-type planet with an envelope (i.e. atmosphere) 
composed of ~85% water by mass mixed into H and He atop a rock core 
with about 0.2 x solar bulk ice to rock ratio. 

Given the vast amount of exoplanets and exoplanet candidates that are 
being discovered at present, we chose to simulate a hypothetical planet 
that has (1) the size, mass, and orbital characteristics of GJ 1214b; (2) 
a greenhouse-type atmosphere (i.e., an atmosphere whose heating/cooling 
properties are dominated by a greenhouse gas without scatterers, akin to a 
dust-free Mars); (3) the gas properties of H2O (but without any sort of surface 
ocean or hydrologic cycle); and (4) a surface, i.e. some density discontinu- 
ity between a gaseous atmosphere and solid or liquid surface that provides 
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a source of friction and interacts thermally with the atmosphere. We in- 
vestigate a range of global mean surface pressures (10 4 -10 7 Pa) and surface 
albedos (0.1-0.7), since these are not yet constrained by observations. 

As shown by our literature review above, exoplanet G CMs are ab u ndant 
We use a state of the art dynamical core comp ared with 



Joshi et al. 



As previously mentioned, 



Showman et al. 



(119971 ) 



(120091 ) use the same dynamical core 



that we do, but use di fferent parameterizations and more importantly lack a 



surface. iMenoul ( 120121 ) has modeled GJ 1214b with a H2O atmosphere, but in 
the mini-Neptune regime since it does not have a solid or liquid surface. This 
is simply a modeling choice; we compare the diffe r ences between our model 



and theirs in Section 



Merlis and Schneider! (120101 ) perform idealized 



simulations, but only for rapidly vs. slowly rotating cases (their rapidly 
rotating case will be useful to compare to our model). 

Heng and Vogt (2011) have already performed paramete r sweeps of super 



Earth planets, but they use Newtonian r elaxation to the 



Held and Suarez 



Held and Suarezl (119941 ) has been 



( 119941 ) radiative equilibrium state (Eq. CQ. 
and continues to be a widely-used, robust, workhorse for terrestrial GCMs. 
Although it has tunable parameters for other planets (e.g., equator to pole 
temperature difference), its functional form is deliberately set so that it im- 
plicitly mimics the features of Earth's atmosphere. We use the same tech- 
nique of Newtonian relaxation, but our expression for the radiative equilib- 
rium temperature is an actual solution to the RT equations (and includes 
convection). Granted, the solution is not completely state-of-the-art in that 
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it is a gray model, but simple models are not without value (more complex 
does not always mean more accurate). Furthermore, no parameter sweep like 
ours has been performed using any of the models with sophisticated RT. 

Another important difference from the Heng and Vogt (2011) study is 
that their default case has a rotation period of 37 Earth-days, whi l e ours 



Merlis and Schneiderl (120101 ) 



assumes the GJ 1214b rate of 1.58 Earth-days, 
have noted the importance in rotation rate, to the extent that the circulation 
is in entirely different regimes for 1 Earth-day periods and 365 Earth-day pe- 
riods. We speculate in Se ction 14.21 that the rotation period of 16 Earth-days 



in the 



Joshietal 



( 119971 ) GCM is also in the slow regime. Heng and Vogt 
(2011) themselves even state that rotation rate is one of the most impor- 
tant parameters that affect the atmospheric circulation. Thus, as Heng and 
Vogt (2011) performed a parameter sweep in the slowly rotating regime, we 
perform a parameter sweep in the rapidly rotating regime. 

In Section [2] we describe the 3-D GCM dynamical core (based on the MIT 
GCM) and external heating/cooling scheme of Newtonian relaxation to a 
prescribed radiative-convective equilibrium temperature profile. In Section [3] 
we present the results for temperature, winds, mass streamfunction, and 
eddy transport, focusing on the 10 5 Pa global mean surface pressure and 
0.4 albedo case. We also perform a comparison of temperature, winds, and 
surface pressure variation as global mean surface pressure and albedo are 
varied. In Section 14.1 



Showman et al. 



we compare our results to the shallow water model of 



(I2009h and other 3-D GCMs in Section EES We also present 
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model transmission spectra in Section 14.31 



2. General circulation model configuration 



2.1. Dynamical core 



Zalucha et al. 



(120101 ) adapted the MIT GCM to Mar s and used it to ana- 



l yze t 



le effect of topography on the Martian Hadley cells. 



Zalucha and Gulbis 



( 120121 ) also adapted the MIT GCM to Pluto and compared the results to ob- 



served stellar occultation light curves. Other groups have also adapted the 



MIT GCM dynam i cal co re to hot Jupiters and hot Neptunes (IShowman et al. 



2009 



Lewis et al. 



2010|). Thus the MIT GCM has been successfully applied 
to other planets beyond its already wide usage as an Earth atmosphere and 
ocean model. 

The dynamical core of the MIT GCM solves t he primitive equations using 



Marshall et al. 



19971 ) on a 



a finite volume method on an Arakawa C-grid 
sphere. The model atmosphere is hydrostatic (see the Appendix for justifi- 
cation of this assumption), compressible, and behaves according to the ideal 
gas law. A thermodynamic equation is included that contains a user-specified 
external heating term (see Section |2~2|) . Convective adjustment is performed 
within the atmosphere to prevent superadiabatic temperature profiles. We 
assume that the gravitational ac celeration is constan t with height. In the 



horizontal, a cubed-sphere grid (lAdcroft et al. 



2004 ) with 32 x 32 points 



per cube face is used, equivalent to a grid spacing of 2.8° or 834 km at the 
equator. 
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The Rossby radius of deformation in the midlatitudes is given by 



where N ~ 0.01 s" 1 is the Brunt-Vaisala frequency, h ~ 10 km is the order- 
of-magnitude scale height, and f ~ 10~ 4 s _1 is the midlatitude Coriolis 
parameter. Then, L# is of order 1000 km. The Rhines length is given by 



where is U is a characteristic velocity (i.e., the square root of the globally 
averaged total kinetic energy), taken to be 100 m s _1 , and (3 is the meridional 
derivative of the Coriolis parameter, taken to be 5 x 10~ 12 (m s) _1 . The 
Rhines length is then approximately 4000 km. Thus, the horizontal model 
grid resolves both the midlatitude Rossy radius of deformation and Rhines 
length. 

Compared to a cylindrical projection grid (i.e., a latitude/longitude grid), 
this type of horizontal grid eliminates singularities at the poles that force 
meridional winds to zero and removes the requirement for Fourier filtering 
in the high latitudes (in order to maintain a practical timestep). Note that 
in contrast to the finite volume method employed by the MIT GCM, other 
exoplanet GCMs (e.g. Rauscher and Menou 2012; Heng and Vogt 2011; 
Cho et al., 2003) use the spectral method, which allows for a cleaner control 
of numerical noise. The MIT GCM configuration has neither viscosity nor 




(3) 
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vertic al diffusion, bu t uses a Shapiro S2 filter to remove purely numerical 



noise ( jShapiro . 



19701 ) . An 8th order filter is used for the momentum. The 
order of the physical space filter is 4. The timescale of the filter is set to 
model time step. Heng, Menou, and Phillipps (2011) have demonstrated 
that numerical horizontal dissipation renders the prediction of temperatures 
to be uncertain at the 10% level and several tens of percent in wind speeds 
in hot Jupiter models; however this was not the case in their Earth models. 
Comparison of our derived wind speeds to that of other studies may suffer 
from this uncertainty. 



The vertical grid is defined in terms of an rj coordinate ( lAdcroft et al. 



2004 ) based on atmospheric pressure and is scaled according to the pre- 
scribed global mean surface pressure, rj coordinates have the advantage of 
being based in atmospheric pressure (which reduce the governing equations 
to a simpler form than height coordinates) and do not suffer large pressure 
gradient errors near the surface when the topography is non-zero (as with 
a terrain-following or a coordinate). In general, surfaces of constant rj will 
intersect the surface, making computation more difficult, though for our pur- 
poses we assume flat topography and this is not a concern. The vertical grid 
has 30 points that vary logarithmically in pressure, such that there are more 
points closer to the surface. Throughout most of the model domain, there 
are about 5-10 points per scale height, except near the bottom where the 
resolution is finer and the top where the resolution is coarser. The top level 
of the model is located at rj = 2 x 10 -5 ; however, the top three points of 
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the model contain Rayleigh drag (discussed below) and are not considered 
physically realistic output. The top most level of the model with physically 
realistic output is located at r\ — 10~ 3 . 

The horizontal momentum equation is given by: 



where D is the material derivative, t is time, is the horizontal component 
of velocity, f — 2Q sin <ft is the Coriolis parameter (where Q is the plane- 
tary rotation rate), $ is the geopotential height, indicates that only the 
horizontal component of the gradient should be considered, and is the 
frictional drag. Here fricti on arises from d r ag on the atmosphere by the 
surface, given by (following iHeld and Suarezl . Il994f l 



where kf is the wind damping rate, p Q is the local surface pressure (a function 
of latitude, longitude, and time), and p b = 0.7p o is the pressure at the top 
altitude of the boundary layer, defined as the part of the atmosphere where 
the drag is nonzero. Eq. |5] represents a drag law (linearly dependent on the 
horizontal velocity) that decreases with height, reaches zero at the top of the 
boundary layer, and is zero at all levels above. It reflects the diminishing 
dynamical influence of the surface the altitude increases and the fact that 
the boundary layer tends to follow the height of the surface terrain (i.e., the 



-fk x u h - V h $ + F h , 



(4) 



Dt 
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top of the boundary layer wi ll 



gener ally be higher where the surface elevation 



Held and Suarezl (119941 ). which is a benchmark GCM for Earth, 

es of 0.2 to 10 



is higher). 

uses kf = 1 Earth-day -1 . Ma rs models use frict ional time sea 



sols 



Joshi et al 



1995|). 1 sol (ILewis et al 



19971 ). and ~20 sols flCaballero et al. 



1996), 0.25-0.7 sol flNavvelt et al 



2008). These relationships are deter- 



mined empirically through modeling of near-surface atmospheric properties. 
Earth and Mars are the only bodies for which empirical values exist, and 
in the lack of any data for exoplanets, we take them to be representative of 
terrestrial-like exoplanets and use kf = 1 day -1 . 



Note that Eq. [5] does not represen t a true 
mimics one using Rayleigh friction. See lGarratt 



jound ary layer scheme, but 



19941 ) for a review of bound- 



ary layer science for Earth. Furthermore note that the method of dealing 
with the momentum at the bottom boundary of hot Jupiter models has also 
not yet reached a consensus. Heng, Frierson, and Phillipps (2011) use a 
free-slip condition, i.e., no friction at the bottom boundary and no condition 



that the velocity para l 



el to the surface should be zero. On the other hand, 



Rauscher and Menoul (120121 ) use Rayleigh friction at the bottom boundary 
to mimic the effects of magnetic drag. 

Rayleigh friction is also used at the model top as a parameterization for 
small-sc ale wa v es th at break (i.e. deposit momentum) in the upper atmo- 



sphere (JLeovyj, 



19641 ). Rayleigh friction is a drag force linearly dependent 
on velocity, similar to the surface friction scheme above, but the force in- 
creases with height in the uppermost layers of the model. Rayleigh friction 
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schemes are tuned to atmospheric measurements where observations of den- 
sity, temperature, and/or wind exists (i.e. Earth, Mars, and Venus), but for 
exoplanets no information regarding small-scale waves is available. Thus, in 
the top three levels we use typical Earth drag coefficients, from top to bottom 
these are 1, 3, and 9 Earth-day -1 . 

Table [T] shows the GCM parameters. We have assumed that the gaseous 
properties are that of a 100% water vapor atmosphere. This choice is some- 
what arbitrary, as we could have chosen another greenhouse-type gas such as 
C0 2 . Nonetheless, we do not expect the fundamental conclusions to be dras- 
tically different. Although we are not explicitly modelin g GJ 1214b, we have 
chosen to use the observed parameters of this planet from 



Charbonneau et al. 



( 120091 ). Along with the assumption of a solid su rface, this combina t ion o f 



Rogers and Seagerl (120101 ) 



parameters represents a hypothetical planet, since 
showed that GJ 1214b cannot have a solid surface if its atmosphere is com 
po sed of water v apor 



Heller et al 



(120111 ) have shown that planets in close orbits around low- 
mass stars will lose their obliquities on very short timescales. Provided that 
putative other planets in the system do not disturb the planet significantly, 
it will not have a significant obliquity, even it was initially close to 90° (R. 
Heller, private communication). Thus, we assume zero obliquity, which by 
definition means the plane t will be at perpetual equinox and have no seasons. 



Heng and Kopparla 



( 120121 ) have shown that the characteristic time scales for 



spin synchronization and circularization of the orbit of GJ 1214b (whose 
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orbital parameters we use) is probably much less than the age of the star. 
Thus, we assume a circular orbit. 

2.2. Newtonian relaxation scheme 

In our GCM the planet is assumed to be tidally locked into a synchronous 
orbit with its parent star. Like the Mars version of our GCM, we use the 
method of Newtonian relaxation to prescribe the external heating/cooling 
term in the model, given by 



where T is the air temperature in the GCM, c p is the specific heat at constant 
pressure, V is volume, T eg c is the radiative-convective equilibrium air tem- 
perature, and krp is the radiative relaxation rate, kx is a parameter that is 
calculated based on the e-folding time of a perturbation to the vertical tem- 
perat ure profile (the method is described in detail in Eq. 6 of lShowman et all 
20081 ). We note that several studies (Heng and Vogt, 2011; Showman and 
Polvani, 2011; Langton and Laughlin, 2007; Thrastarson and Cho, 2011) have 
pointed out the importance of radiative relaxation rate on the atmospheric 
circulation, but assume constant fcr here for simplicity. 

In deriving the expression for T e(?)C , we first start with the radiative equi- 
librium temperature T eq (and its value at the surface T eg )S ) that does not 
contain the effects of convection; we refer to it as the pure radiative equi- 
librium temperature. For an atmosphere that is gray at long wavelengths 





24 



and transparent at short wavelengths, the RT equation in the Eddington 
approximation has the simultaneous analytical solution for the atmosphere 



aT e 4 9 = (l-A)Q(0.5 + 0.75r) (7) 

and surface 

aT^ s = (1 - A) Q (1 + 0.75r s ) , (8) 

where a is the Stefan-Boltzmann constant, A is the surface albedo, r is the 
optical depth, r s is the optical depth at the surface, and Q is the stellar flux 
received by the planet at a particular latitude phi, hour angle h (related to 
longitude), and stellar declination 5 (related to season). Q is given by, 

Q = Qo cos (sin sin 5 + cos cos 5 cos h) , (9) 

where Q Q is the stellar flux received by the planet (equivalent to the so- 
lar constant defined for Earth). This radiative equilibrium model predicts 
the surface temperature T eg s by flux balance with the outgoing surface flux 
(crT e 4 ), the downwelling long wave flux from the atmosphere, and the down- 
welling short wave flux from the star. As such, the only assumption about 
the properties of the surface is that it have a density discontinuity with the 
atmosphere (e.g., is a liquid or solid). The surface is infinitesimally thin and 
has zero heat storage capacity (i.e. no thermal lag) in this radiative equilib- 
rium model (although we could account for an emissivity not equal to unity 



25 



by multiplying the left hand side of Eq. El by the emissivity) . Note that the 
surface temperature of the GCM output may be prognostically determined 
using a slab method or other method of the user's choosing. 

The optical depth is linearly dependent on pressure as r = Tn„(p/pn n), 



20101 ) is 



where t d o (whose numerical value is taken from iMer lis and Schneider! , 
the optical depth at the reference pressure p OQ = 10 5 Pa. Thus, as the surface 
pressure is varied, so will r at the surface, indicating a stronger greenhouse 
effect for more massive atmos pheres. Our gr ay atmosphere specification is 



equivalent to that derived by 



Mihalasl (119781 ) in the limit that the interior 



radiation term (i.e., infrared radiation emitted by the surface) dominates 
the shortwave stellar radiation directly absorbed by the atmosphere. Note 
that at higher pressures (> 10 6 Pa), the assumption that the atmosphere 
is transparent to short wavelengths may not be physically valid, and the 
lower levels of these thicker atmospheres may be optically thick at short 
wavelengths. 

Note that, within the framework of our model, it is possible to rewrite r 



as r = TnJv/Van) n , where n = 2 



ity (ILorenz and McKay , 
tiona 



20031 ) . 



represents a source o 



collisi on-induced opac- 



Merlis and Schneider! (120101 ) incl ude an addi 



Heng et al. 



term that is proportional to the fourth power of pressure. 
( 120121 ) have developed a more complicated semi-analytical formalism for cal- 
culating the temperature-pressure profiles in hot Jupiter atmospheres with 
a correction factor for collision-induced opacity. Collision-induced opacity is 



an avenue of future improvement to our model. 
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In the mathematical limit r — > t s in Eq. there is a temperature dis- 
continuity between the pure radiative equilibrium temperature of the surface 
and the pure radiative equilibrium temperature of the atmosphere just above 
the surface. To mitigate this unphysical property of th e solution, convective 



adjustment is performed (see 



Lorenz and McKay 



2003 



Zalucha et al. 



2010 . 



for details of this scheme) on the pure radiative equilibrium atmospheric and 
surface temperatures. In radiative- convective equilibrium, the temperature 
profile within the near-surface convective layer follows an adiabat and is con- 
tinuous with the surface temperature. Explicitly, the radiative-convective 
equilibrium temperature is 



eq.c 



[l -A) Q (0.5 + 0.75r t ) 



AR/c p 



(10) 



for r > Tt, where r t is the r at the top of the convective layer. For r < r t , T eQ}C 
is given by Eq. [7J r t is solved by numerically balancing the upward surface 
flux in pure radiative equilibrium, the integrated upward atmospheric flux 
for r > r t in pure radiative equilibrium, the upward surface flux in radiative- 
convective equilibrium, and the integrated upward atmospheric flux in the 
convective layer in radiative-convective equilibrium. Specifically, 



° T Ls exp [- (t s -T t )/fi]+ / crT e 4 (r*) exp [- (r, - r t ) ///] dn 



(11) 



° T LA T s) exp [- (r s - T t ) In] + / aT 4 e{ ' (r*) exp [- (r* - r t ) ///] dn 
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where /x = 2/3 is the average cosine of the emission angle and r* is a variable 
of integration. T e(JjC is now quantitatively known for use in Eq. |6j 

Under the condition of surface- atmosphere mass exchange, the atmo- 
sphere deposits ice on the surface (if the atmospheric temperature drops 
below the condensation temperature) or the surface sublimates ice into the 
atmosphere (if the surface temperature is greater than the condensation tem- 
perature). In both cases the atmosphere remains at the condensation tem- 
perature. For arguments of simplicity, we do not include surface-atmosphere 
mass exchange (but this is an avenue of future work). We do account 
for the temperature exchange by prescribing that if the radiative equilib- 
rium or GCM-predicted temperature drops below the condensation tem- 
perature (such as on the nightside of the planet where the insolation is 
zero), the temperature is instantaneously reset to the condensation tem- 
perature. Otherwise, the temperature could drop to absolute zero, which 



is unphysical. T his convention has 



Deen demonstrated to be valid in the 



Mars community (jHaberle et al 



1993 



Takahashi et al 



2003 



Zalucha et al 



2010), except near the poles where surface-atmosphere mass exchange oc- 
curs. This step implicitly supplies the latent heat that would have been 
added/subtracted, thus conserving energy. The condensation temperature is 
given by the Clausius-Clapeyron relation 



T, 



cond 



1 R, 
In 

Ti L 



p_ 

J*. 



-i 



(12) 
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where T\ and p\ are the vaporization reference temperature and pressure, 
respectively (see Table [[]) and L is the specific latent heat of vaporization. 
Note that L has a dependence on temperature, which we do not take into 
account (since we are not aware of any measurements taken at the high 
temperatures encountered on hot exoplanets). 

3. Simulation results 

We performed simulations starting from rest, initializing with a globally 
constant temperature and surface pressure. Figure [[] show s the globally aver- 



aged kinetic energy per unit mass for selected cases (see e.g. 



1969 



Leovy and Mintzl . 



for a demonstration of this method). After 120 Earth-days, the model 
reached a quasi-steady state. The results presented for the remainder of this 
paper are the time average of the subsequent 400 Earth-days. A 120 Earth- 
day spin up ti me is short compared to the 200 Earth-day spin up time of 



the benchmark 



Held and Suarezl (119941 ) Earth model. We suspect that this 



is because our radiative relaxation time is shorter than their 40 Earth-day 
value. 

The globally averaged surface pressure, which is a representation of the 
atmosphere's mass, is an important but unconstrained parameter for super- 
Earths. Since our atmosphere is assumed to be entirely made up of the 
radiatively active constituent, the global mean surface pressure is directly tied 
to the strength of the greenhouse effect. Higher global mean surface pressure 
means more IR absorption, and thus higher temperatures. We have chosen 
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to perform simulations using the following global mean surface pressures: 
10 4 , 10 5 , 10 6 , and 10 7 Pa. For reference, the global mean surface pressures 
of Venus, Earth, Mars, and Titan are approximately 90 x 10 5 Pa, 10 5 Pa, 
0.006 x 10 5 Pa, and 1.5 x 10 5 Pa, respectively; however, note that none of 
these atmospheres are purely greenhouse-gas dominated as we have assumed 
in our study. 

The geometric albedo of a planet is also an important parameter for plan- 
etary atmospheres because it controls the amount of energy that is absorbed 
by the planet's cross-section, and thus the amount of energy absorbed by the 
surface and re-emitted into the atmosphere. Albedo variations (e.g. due to 
inhomogeneous surface distributions of oceans, rocks, and/or ices) over the 
planet's surface can also be important. For simplicity, we perform simula- 
tions with globally constant surface albedos of 0.1, 0.4, and 0.7. Note that 
the globally averaged surface albedos of Venus, Earth, Mars, and Titan are 
approximately 0.77, 0.3, 0.15, and 0.22, respectively. 

3.1. Results of the 10 5 Pa surface pressure, 0.4 surface albedo case 

We examine in detail the results of the 10 5 Pa global mean surface pressure 
and 0.4 surface albedo case. This case is near the center of our parameter 
sweep and is thus a prudent starting point for our discussion. 

Figure [2] shows the zonally averaged potential temperature. This plot 
and the subsequent zonal mean plots are plotted with a level = p/p s as the 
vertical coordinate (i.e., the plots have been normalized to the global mean 
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surface pressure). At a levels greater than 0.02 (pressure 2 x 10 3 Pa), at a 
fixed level, there is a dip in potential temperature in the tropics. At lower a 
levels, there is a rise in potential temperature in the tropics. 

Figure [3] shows temperature profiles from the equator and 60° latitude 
(midlatitude) at four local times. All of the temperature profiles contain a 
temperature inversion (temperature increasing with decreasing pressure) near 
the surface, and the midlatitude profiles contain another inversion peaking 
at about the 0.01 a level (pressure 10 3 Pa). The near-surface inversions are 
also present in some of the other global mean surface pressure and surface 
albedo cases. The precise cause of this inversion is unknown. The underlying 
radiative-convective temperature monotonically decreases with decreasing 
pressure (Eq. [7]) because there are no assumed shortwave absorbers and the 
optical depth decreases linearly with pressure. Moreover, on the nightside, 
where the insolation is zero (and therefore the radiative-convective tempera- 
ture would be unphysically at absolute zero), the column radiative-convective 
temperature is set to the condensation temperature. Since the radiative re- 
laxation timescale is assumed constant with height, it is puzzling that the 
surface temperature would be pulled more strongly towards colder radiative- 
convective equilibrium temperature and yet at a slightly higher altitude, less 
so. 

Temperature inversions may be caused dynamically be two physical pro- 
cesses: downwelling, which leads to adiabatic compression/heating, and dif- 
ferential horizontal heat transport with height. We have searched through 
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all of our cases and cannot find a direct correlation between the presence 
of an inversion and downwelling. We have not analyzed the local horizontal 
heat transport, as this is beyond the scope of this paper. We speculate that 
the cause of the temperature inversion is some complex combination of these 
two physical processes. 

Comparing pressure-temperature (p-T) profiles at the same local times 
indicates that at some altitudes the temperature is higher in the midlatitudes, 
which would not be expected through astronomical considerations alone be- 
cause the stellar flux decreases as a function of latitude. The reversal occurs 
at pressures lower than about the 0.01 a level (10 3 Pa) and within a few a 
level tenths (10 4 Pa) of the surface. There is also a longitudinal asymmetry 
between the morning and evening terminators. At the equator, the evening 
temperature profile is greater than or equal to the morning temperature at 
all levels, while in the midlatitudes, the opposite is the case. Note that the 
temperature profiles are always greater than the condensation temperature 
of H2O, because we have explicitly forbidden it in the model equations and 
parameterizations. 

Figure H] shows the wind structure of the free atmosphere (above the 
frictional boundary layer) at a level 0.068 (pressure 6800 Pa). While quanti- 
tative differences exist at the levels outside the boundary layer (pressures less 
than 0.7 x 10 5 Pa), they all exhibit a stationary wavenumber-1 longitudinal 
structure in the vertical velocity field. The circulation is divided into two 
latitudinal regimes. First, an equatorial region between ±30° with a high 
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speed (> 300 m s _1 ) westerly jeto, sinking motion on the western (i.e., morn- 
ing) hemisphere, and rising motion in the eastern (i.e., evening) hemisphere. 
Second, poleward of ±35° latitude, the flow is northward over the pole on 
the terminators and dayside and eastward around the nightside of the planet. 
Moreover, the sign of the vertical motions switch such that descent is on the 
evening hemisphere and ascent is on the morning hemisphere. 

Figure |5] shows the zonally averaged zonal mean wind. There is a slight 
dependence on altitude, but westerly winds occur in the tropics between 
30° to 40° in each hemisphere. Easterly winds occur at higher latitudes. The 
maximum westerly wind strength is achieved at around a level 0.1 (pressure 
10 4 Pa), while the maximum easterly wind strength occurs at around a level 
0.002 (pressure 200 Pa). Figure [6] shows the latitudinal transport of westerly 



eddy momentum (per unit mass), u'v', where u is the zonal velocity, v is 
the meridional velocity, primes denote deviations from the zonal means, and 
the overbar denotes a zonal average. This figure shows that in the Northern 



hemisphere, u'v ' < 0, indicating southward transport of westerly momentum, 



while in the Southern hemisphere u'v' > 0, indicating northward transport. 
Thus, westerly eddy momentum is converging at the equator. The magnitude 
of the transport increases with height at all latitudes, and is strongest at 
approximately ±45°. 

The two latitudinal regimes seen in the zonal wind fields may also be seen 



1 For reference, the sound speed in an ideal gas of water vapor is 495 m s 1 at 400 K 
and 858 m s" 1 at 1200 K 



33 



in a plot of zonally averaged Euler-mean mass stream function (Fig. [7]). In 
each hemisphere, there is a cell that rises at the equator, flows poleward aloft, 
and sinks at ~ 35° latitude. Poleward of this latitude, there is a cell that 



circulates in the opposite sense and is s tronger than the equat o rial fl ows. The 



one Earth-day rotation period case of 



Merlis and Schneider! (120101 ) shows a 



three celled structure. Moreover, their maximum stream function magnitude 
is two orders of magnitude lower than ours. This is consistent with our 
higher temperatures and zonal and meridional wind speeds. The atmosphere 
is more energetic and must transport heat away from the substellar point 
more vigorously. 

In the boundary layer at a level 0.9377 (pressure 93770 Pa), the circula- 
tion has a different structure than in the free atmosphere. The midlatitudes 
are characterized by weak downwelling and a horizontal wavenumber-1 pat- 
tern. At the equator, an intense area of downwelling and easterly flow is 
present at the substellar point (12 hours LT), while there is an area of in- 
tense upwelling at 15 hours LT. Still farther to the east the upward motion 
switches back to downward. 

3.2. Results of differing surface pressure and surface albedo 

The global mean surface pressure (i.e. the atmosphere's total mass) and 
the surface albedo are not well constrained by observations; however, these 
are fundamental parameters that strongly modulate atmospheric processes. 
In order to compare simulations with different global mean surface pressures, 
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subsequent quantities are plotted on the same a level . If we instead chose a 
constant pressure level, we might be comparing a level that is in the boundary 
layer of a low surface pressure case with a level that is near the model top in 
a high pressure case, where different physical processes are acting. 

Figure [8] shows dayside, equatorial temperature profiles as a function 
of p s and A for both the GCM simulated temperatures and the radiative- 
convective equilibrium temperatures (see Section 1272]) . As p s increases, tem- 
perature increases. This behavior occurs simply because we have assumed 
that the atmosphere is entirely composed of the IR absorber, so as atmo- 
spheric mass increases, so does the amount of the absorber. The atmosphere 
subsequently absorbs more of the energy being reradiated by the surface. 
Similarly, temperature decreases as A increases because more of the incoming 
stellar flux is reflected to space. In general, the GCM simulated temperatures 
are cooler than the radiative-convective equilibrium temperatures (for a given 
p s and A) in the lower atmosphere, but warmer in the upper atmosphere. 
The GCM temperatures are more isothermal. By definition, the radiative 
equilibrium temperature is the temperature of a motionless atmosphere, so 
any difference is indicative of circulation. In Fig. [HJ adiabatic heating due 
to downward motion is either occurring on the precisely at LT=12 hours 
and the equator, or at another location with subsequent lateral transport to 
LT=12 hours at the equator. 

These may be lateral heat transport by horizontal winds and/or adiabatic 
heating and cooling due to downward and upward motions, respectively. 
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Figure M shows GCM temperature in the free atmosphere (outside of the 
frictional boundary layer) as a function of p s with A held constant. As p s 
increases, the temperature increases at all locations, for reasons described 
above. In all cases, the equatorial temperature maximum is shifted east- 
ward of the substellar point, as is also t he case in hot Jupiter models (e.g., 



Showman et al. 



2008 



Perna et al. 



20121 ) . The maximum temperature exists 



offset from the equator at about 30° latitude and 12 hours LT. Thus, on the 
dayside as we move from equator to pole, the temperature first increases until 
this location is reached, then decreases towards the poles. At high latitudes, 
the temperature on a given latitude circle is maximum near 12 hours LT, but 
the exact longitude of the maximum is located on either side of this merid- 
ian depending on the direction of the zonal flow. A temperature minimum 
occurs at hours LT and 50° latitude. On the nightside, the temperature 
decreases monotonically from equator to pole. 

Figure [TO], along with panel (b) of Fig. [9] (in the order Fig. [TUb . Fig. [9b, 
Fig. dOb), show the effect of increasing A on temperature in the free atmo- 
sphere (p s held constant). Again, the effect on the horizontal structure is 
simply to decrease the temperature everywhere as A increases. 

Figure (TT] shows the zonal component of wind in the free atmosphere as 
a function of p s . As p s increases, the speed of the equatorial jet decreases 
between 9 and 15 hours LT. The latitudinal extent of the equatorial jet 
increases on the nightside. At the two greatest initial surface pressure cases 
(10 6 and 10 7 Pa), an increasingly westerly component of flow develops at 
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middle to high latitudes between 11 and 18 hours LT. The nightside, easterly 
flow in the midlatitudes also increases in speed as p s increases. Figure [T2~| 
along with panel (b) of Fig. [TT] (in the order Fig.lT2"a. Fig.fTTb. Fig.[T2"b). show 
the effect of increasing surface albedo on zonal wind in the free atmosphere. 
The zonal wind speed in the equatorial and midlatitude jets decreases with 
increasing A. Also, the latitudinal width of the equatorial jet decreases as A 
increases. 

Figures [13] and [14] show the effect of increasing p s and A, respectively, 
on meridional flow in the same order as Figs. [TT] and [T2J All cases have 
areas of southerly flow at southern longitudes from 12-24 hours LT and 
northern longitudes from 0-12 hours LT and northerly flow elsewhere. In the 
cylindrical projection format of our plots, this behavior indicates that there 
is a component of the flow over the poles from the evening side (i.e., western 
hemisphere) to the morning side (i.e., eastern hemisphere). As p s increases, 
the meridional wind speed increases; as A increases the speed decreases. 
This behavior implies a link between the amount of thermal energy and the 
strength of the winds. 

Figure [15] shows the time-integrated surface pressure as a function of 
latitude and longitude and p s . Note the distinction between the global mean 
surface pressure (p s ) and the latitudinally and longitudinally varying surface 
pressure — the latter is equal to p s in the global sum. To show each case on 
the same scale, the latitudinally and longitudinally varying surface pressure 
has been normalized by p s for that case, so what is actually shown is the 
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surface pressure anomaly. The surface pressure is higher at the poles than 
at the equator, which is opposite of terrestrial-type Solar System bodies 
(Earth, Mars, Venus, and Titan). As p s increases, the fractional difference 
between pole and equator decreases. In the p s = 10 4 and 10 5 Pa cases, the 
surface pressure variation is more longitudinally symmetric, with an entirely 
positive anomaly poleward of ~ 40 — 60° and an entirely negative anomaly 
equatorward of this latitude circle. In the p s = 10 6 and 10 7 Pa cases, the 
surface pressure anomaly is entirely positive poleward of this latitude circle, 
but equatorward there are both positive and negative anomalies for a given 
latitude circle. The minimum surface pressure anomaly moves poleward and 
westward as p s is increased; the maximum surface pressure anomaly moves 
mainly equatorward. 

Figure [161 along with panel (b) of Fig. [15] (in the order Fig. [T6b . Fig. [15b, 
Fig. fl~6b). show the effect of increasing surface albedo on the surface pressure 
anomaly field. The minimum negative anomaly deepens and moves poleward 
and westward, but otherwise the effect of changing surface albedo is small 
on the surface pressure field. 

Within the boundary layer (here for the case of a level 0.9377), as in 
the free atmosphere, temperature increases as p s increases. The temperature 
maximum is located on the equator and at the substellar point (12 hours 
LT), although the contours are elongated in the eastward direction. The 
temperature minimum is also located on the equator at hours LT, except in 
the p s = 10 7 Pa case where it is at ~ 40° latitude. As in the free atmosphere, 
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temperature increases as A decreases. 

In the boundary layer, a localized area of westerly wind exists between 
about 3 and 15 hours LT and 40° latitude. The maximum speed in this 
area decreases with increasing p s and A. Poleward of 40° latitude there is an 
easterly jet that is most prominent in the low p s and A cases. As both p s and 
A increase, the jet strength decreases and the maximum speed shifts locations 
to the equator at 18 hours LT. The meridional wind field is complicated, but 
is generally a mirror image of itself reflected across the equator. A clear trend 
does not exist as p s and A are increased, except that in the p s = 10 7 case 
there is an abrupt amplification of wind speeds. 

4. Discussion 



4-1. \Showman and Polvani (201 A) analytic shallow water model 

We pay special attention to the lShowman and Polvanil (120111 ) shallow wa- 



ter analytic models, as they have already laid the analyt ical groundwork that 



will he lp explain our model results in si mple r term s. The 



(120111 ) model is a generalization of the 



Gill fll980h and 



Showman and Polvani 



Matsund (119661 ) mod- 



els to cases where the radiative and frictional timescales are not equal, as is 
the case for tidally locked, synchronously orbiting, strongly irradiated, short- 
period exoplanets. This type of model leads to a solution where the day-night 
thermal forcing generates standing, planetary-scale Rossby and Kelvin waves, 
which interact with the mean flow. The Kelvin waves span the equator, 
and the Rossby waves lie at the poleward boundaries of the Kelvin waves. 
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The combination of westerly propagating Kelvin waves; easterly propagating 
Rossby waves; and the pressure gradient, Coriolis, advective, and drag forces, 
the velocity solutions resemble eastward pointing "chevrons" at the equator. 
These velocity tilts transport westerly momentum from high latitudes to the 
equator, creating the equatorial superrotation. 

We find that the geopotential (i.e. temperature) field of their non-linear 
solution with radiative and drag constants of one day most closely matches 
our temperature pattern qualitatively (see Figs. [9] and [10]); that is, an east- 
ward phase shift of the temperature at the equator and a westward phase shift 
at high latitudes. The temperature contours exhibit northwest-southeast 
tilts in the Nort hern hemisphere and southwes t-northeast tilts in the South- 



ern hemisphere. 



Showman and Polvanil ( 120111 ) describe this as an eastward 



pointing chevron pattern. Furthermore, the absolute temperature maxima 
and minima are located off the equator. This behavior can be seen in the zon- 
ally averaged potential temperature plots (Fig. HJ at a level 0.068 (i.e., the a 
level of Figs. 191 and fTUl) . where, scanning from pole to equator, the potential 
temperature first increases slowly, then drops more sharply in the vicinity 
of the equator. Note that a more appropriate case for comparison with our 
GCM would be a radiative time constant of 10 days and drag constant of 



one day, but 



Showman and Polvanil ( 1201 ll ) do not show this case 



The equatorial Rossby radius of deformation is given by 



1/2 



■iR 



ft 



(13) 



40 



where g is the gravitational acceleration at the surface, d is the shallow- 
water approximation depth (taken to be the depth of the model below the 
Rayleigh spong e layer, 10 3 km) and = 2Q co s (f)/a, where a is the planetary 



radius. In the 



Showman and Polvanil (120111 ) model, the latitudinal extent 



of the equatorial jet is ±Lr. In our GCM, Lr = 6550 km or 22°, which is 
smaller than the jet width of ~ 30-45° (see Fig. E]). Note that the length 
scale Ln is well resolved by our h orizontal model grid. 



One key difference between the 



Showman and Polvani (2011) shallow wa- 



ter m odel and our GCM (and hot Jupiter models, see 



2011 



Showman and Polvanil . 



and references therein) is that in our GCM, the winds are westerly at 



all longitudes at the e q uator , not divergent from near the substellar point. 



Showman and Polvanil ( 120111 ) state that the shallow water (i.e., two layer) 



model does not capture the full multilayer structure of a 3D atmosphere. 



hey do however note that during the spinup phase of the 



Showman et al. 



( 120091 ) hot Jupiter GCM, the equatorial wind field is divergent like the shal- 
low water model, before transitioning to the fully spun up case where the 
zonal winds are westerly at all equatorial longitudes. We find precisely the 
same behavior. At 30 days of integration, our model shows the divergent 
zonal wind pattern — easterly winds from to 12 hours local time and west- 
erly winds from 12 to 24 hours local time. At 60 days of integration, the 
region of easterly winds at the equator has shrunk to 6 to 15 hours local 
time, and by 90 days (which from Fig. [1] is now entering the fully spun up 
phase), the equatorial winds are entirely westerly. 
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4-2. Comparison with previous exoplanet GCMs 

Here we focus primarily on a comparison with super-Eart h/GJ 1214b 



mode 



1997 



s that most closely resemble our mode l configuration (IJoshi et al. 



Merlis and Schneider 



2010 



Menou 



20121 ) and parameter sweep stud- 



ies (Heng, Menou, and Phillipps, 2011; Heng and Vogt, 2011). These models 
have been already summarized in Section 11.11 



Our resu 



Joshi et al- 



ts do not compare well with the terrestrial exoplanet GCM of 



(119971 ). if we consider the 10 5 Pa surface pressure cases (their 
Run 1). Their GCM produces one overturning cell in each hemisphere, 
whereas our GCM produces 2. Their temperature maximum is at the equa- 
tor and the 0.8 a level. Their temperature decreases monotonically towards 
the poles; our temperature maxima are offset from the equator. The tem- 
perature magnitudes do not match, although this may be accounted for by 
the fact that we are using a much higher incoming stellar flux. Their zonal 
wind pattern shows westerly flow everywhere, that increases with height and 
decreases from equator to pole. Meanwhile, our GCM shows the distinct 
equatorial westerly jet at the equator and easterly flow at high latitudes. 
Their stationary horizontal eddy momentum flux (equivalent to our Fig. [6]) 
has much more structure, whereas ours are simply negative in the Northern 
Hemisphere and positive in the Southern Hemisphere. We speculate that 
these differences are due to their assumed rotation period of one Titan day 
(about 16 Earth-days) , comp ared with our rotation period of 1.58 Earth- 



days. The I Joshi et al. 



(119971 ) more closely resembles in a qualitative sense 
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the slowly rotating case of 



Merlis and Schneider! ( 120101 ). which has one over 



turning cell in each hemisphere. 



Merlis and Schneider! (120 lOf ). with a rotation 



The quickly rotating case of ] 
period of 1 Earth-day, 10 5 Pa surface pressure, and surface albedo of 0.38 has 
the most similar configuration to our 10 5 Pa surface pressure and 0.4 surface 
albedo case. Their surface temperature pattern, while not being a quan- 



titative match since they assumed te rrestrial levels of insola t ion, 



the chevron pattern of our GCM and 



Showman and Polvanil (120111 ). having 



ooks like 



the maximum offset from the equator and a local minimum on the equa- 
tor. However, their zonal winds at high altitudes are actually easterly at the 
equator (except for weak westerlies just east of substellar point) and westerly 
at high altitudes. Thi s patt ern is completely opposite from our GCM and 



Showman and Polvanil (120 111 ). Furthermore, the 



Merlis and Schneider! (120101 ) 



GCM results have three overturning cells in each hemisphere, compared with 
two in our GCM. We suggest two possible causes for this discrepancy: ei- 
ther the slightly faster rotation rate causes the atmospheric circulation to 
enter a different regime (it is a general rule for Solar System planets that 
faster rotation rates produce more overturning cells) or there is feedback 
with the precipitation scheme in their model. While they do not include the 
radiative effects of clouds, latent heating is present that may be adding an ad- 
ditional heating or cooling source to the model. Furthermore, the winds may 
transport water vapor; thus the circulation and latent heating are strongly 
coupled. 
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The configuration of the GCM of 



Menoul (120121 ) is slightly different than 



our GCM in that they do not include a surface density discontinuity; i.e., their 
model of GJ 1214b is a mini-Neptune class rather than a super-Earth class. 
They assume an entirely H 2 vapor atmosphere with the lowest layer depth 
at 10 6 Pa. Their temperatures are slightly cooler than ours, and do not have 
an inversion in the lowest layer. Their equatorial temperature profiles are 
isothermal to about the 0.05 a level. These differences are almost certainly 
due to the different external heating/cooling schemes (i.e. RT vs. Newtonian 
relaxation as well as the specific dependence of opacity on pressure) between 
the two models. Their zonally averaged winds are westerly everywhere except 
at the deepest levels. They do have an equatorial jet, but at high latitudes 
the winds are still westerly, compared with easterly in our GCM. This is 
possibly a secondary consequence of the different external heating/cooling 
schemes, but may also have to do with the drag at the bottom levels. They 
parameterize magnetic drag while we parameterize frictional drag, but in any 
case their drag time scale is 20 and 40 days in the lowest two model levels, 
while ours decreases from 1 day at the lowest level to at the 0.7 a level. 



Held and Suarezl (119941 ) bench 



Heng, Menou, and Phillipps (2011) use the[ 
mark as their equilibrium temperature specification, modified to include the 
geometric effects of tidal locking and synchronous rotation (Eq. [[]), to simu- 
late a hypothetically tidally locked and synchronously orbiting Earth (surface 
pressure 10 5 Pa). The y use a rotation rate equal to one Earth- year, as in the 



slowly rotating case of 



Merlis and Schneider! (120101 ) . Not unexpectedly, their 
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results do not match ours, having a surface temperature maximum located 
at the substellar point and zonal winds that are divergent from the substellar 
point. 

Heng and Vogt (2011) perform a parameter sweep over many different 
GCM parameters; the one in common with our parameter sweep is the global 
mean surface pressure. Heng and Vogt (2011) u se Newtonian relaxation 



Held and Suarez 



(1994) 



to the equilibrium temperature specification from 
(modified to include the geometric effects of tidal locking and sychronous 
orbit, Eq. [Q. Their rotation rate is 37 Earth-days, compared to 1.58 Earth- 
days in our si mulations. Their r esults for zonally averaged winds are very 
reminiscent of 



JoshietaL 



( 119971 ) Run 1; that is, westerly everywhere, in- 
creasing with height, and decreasing from equator to pole. Heng and Vogt 
(2011) conclude that there is no qualitative difference in the zonal mean zonal 
wind profiles for the cases they tested (10 4 , 10 5 , 3 x 10 5 , and 10 6 Pa). We 
speculate that this is because they are in the slowly rotating regime, while 
we are in the rapidly rotating regime, and it is this point that makes our 
parameter sweep novel. 

We finally note that our GCM shares the entirely westerly equatorial jet 
shown by certain other fully spun up hot Jupiter and super Earth GCMs (e.g. 
Showman and Guillot, 2002; Cooper and Showman, 2005; Showman et al., 
2008, 2009; Heng, Menou, and Phillipps, 2011; Rauscher and Menou, 2012). 
This seems to be a robust feature of synchronously ro tating exoplanets due to 



their subsequent intense day-night heating gradient ( jShowman and Polvanil . 
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201ll ). with the caveat that they must be quickly (with periods of order 1 
Earth-day) rotating. It may also be due to the fact that our model domain 
extends much higher than most Earth or super-Earth models. 



4-3. Transmission spectra 

We compute model atmospheric spectra from our GCM output to deter- 
mine if the atmospheres we have simulated here could be distinguished by 
observations. Specifically, we compute tr ansmission spectra using a l ine-by - 



Madhusudhan and Seagerl (120091 ). 



line radiative transfer code adapted from 
Given the pressure-temperature (P-T) profiles at the limb of the atmosphere 
and the assumed H 2 composition, our spectral model computes a transmis- 
sion spectrum of planet under the assumption of local thermodynamic equi- 
librium and hydrostatic equilibrium. The P-T profile in our spectral model 
comprises of 100 layers uniformly separated in pressure between p s and 1 Pa 
(i.e, 10~ 5 bar), interpolated between the pressure layers of the P-T profiles 
we adopt from the GCM output. In cases where our GCM output does 
not extend to low enough pressures, the atmosphere is assumed isothermal 
with the temperature set equal to the temperature of the uppermost GCM 
level below the Rayleigh friction layer. This is roughly valid for two reasons. 
First, at low pressures the temperature profile is expected to asymptote to 
an isotherm (i.e., skin temperature). Second, the actual temperature may be 
off by ~100 K or less from the correct profile, but transmission spectra for 
exoplanets are presently insensitive to such temperature differences. 
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Figure [T7| shows our model transmission spectra for the different cases we 
consider in this work. Infrared spectra are shown for different global mean 
surface pressures and surface albedos. As shown in Fig. (TTJ the infrared 
transit depth increases with increasing p s over the entire wavelength range 
between 1 and 10 /im. This is simply because with increasing surface pres- 
sure the infrared photosphere (the radius below which effectively all incident 
starlight is blocked from the observer) of the planet is reached at a higher 
altitude, causing a higher effective radius of the planet. Current observa- 
tions would not be able to distinguish variations in surface albedo of our 
hypothetical planet. Observations with existing instruments may however 
be able to distinguish between atmospheres varying in surface pressure by 
over a decade (e.g. between 0.1 and 10 bar, 1 and 100 bar, and 10 and 100 
bar). As shown in Fig. [T7J such differences in surface pressure translate to 
>200 ppm differences in transit depths across the near to mid infrared wave- 
lengths, and hence are detectable with current instrume nts which routinely 



obtain precisions of 100 ppm and less (jBean et al 



201 lb 



The differences between the spectra from models that differ by more than 
two decades in p s are observable with current instruments (i.e. 0.1 and 
10 bar, 1 and 100 bar, etc.). For higher p s , even models with p s = 10 bar 
and p s = 100 bar are distinguishable. The variation in the infrared transit 
depth as a function of albedo is smaller than the variation with pressure. 
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5. Conclusions 



We have simulated, using a GCM with a Newtonian relaxation scheme, a 
hypothetical, tidally locked, synchronously orbiting super-Earth under the 
following assumptions: (1) the size, mass, and orbital characteristics of 
GJ 1214b, (2) a greenhouse-gas dominated atmosphere, (3) the gas prop- 
erties of H2O, and (4) a surface. We have performed a parameter sweep in 
global mean surface pressure (10 4 , 10 5 , 10 6 , and 10 7 Pa) and global mean sur- 
face albedo (0.1, 0.4, and 0.7). Our study differs from the parameter sweep 
of Heng and Vogt (2011) in that we are in the rapidly rotating planet regime 
(diurnal period of 37 Earth-days vs. 1.58 Earth-days), which has been explic- 



itl y shown to be marked 



by 



y different from the slowly rotating regime examined 



Merlis and Schneider! (120101 ) . The other difference in our study from Heng 



and Vogt (2011) is that our relaxation scheme relaxe s to a true analyti c al (al- 



beit gray) solution of the RT equation, ins t ead of the 



specification. While the 



Held and Suarezl (1199 



Held and Suarezl (119941 ) equilibrium temperature is 



the golden standard for benchmark studies of the terrestrial atmosphere, it 
is an intentionally prescribed function that reproduces Earth's atmosphere, 
and therefore may not be general enough for other planets. 

An equatorial, westerly, superrotating jet is a robust feature in our GCM. 
At high latitudes, the flow is easterly. One key point is that the equa- 
torial jet is westerly at all longitudes. This may be reconciled with the 



Showman and Polvanil (2011 ) 



the 



Gilll dl980|) and 



analyt ical models, which is a generalization of 



Matsunol (119661 ) shallow water models, by noting that 



48 



during spinup phase, full 3-D GCMs (e.g. IShowman et al.l . 120091 ) do show di- 
vergent zonal flow at the equator (interpreted as easterly and westerly prop- 
agating Rossby and Kelvin waves, respectively). However, when the steady 
state is reached, the equatorial flow becomes entirely westerly and the anal- 



ogy w ith the shallow water model breaks down, which IShowman and Polvani 



(120111 ) suggest is due to the multi-layered configuration of 3-D GCMs. Our 
GCM results reproduce this behavior precisely. 

Unlike Heng and Vogt (2011), our zonal winds do show a change with 
global mean surface pressure. As global mean surface pressure increases, 
the speed of the equatorial jet decreases between 9 and 15 hours LT. The 
latitudinal extent of the equatorial jet increases on the nightside. For the 
two greatest initial surface pressure cases (10 6 and 10 7 Pa), an increasingly 
westerly component of flow develops at middle to high latitudes between 
11 and 18 hours LT. The nightside, easterly flow in the midlatitudes also 
increases in speed as global mean surface pressure increases. Furthermore, 
the zonal wind speed in the equatorial and midlatitude jets decreases with 
increasing surface albedo. Also, the latitudinal width of the equatorial jet 
decreases as surface albedo increases. 

This study, i.e., a parameter sweep in the rapidly rotating super-Earth 
regime, could be repeated after several model improvements such as replac- 
ing the Newtonian relaxation scheme with directly calculated heating and 
cooling rates, a state-of-the-art (e.g. double-grey, dual-band, or correlated- 
k) radiation scheme, surface/subsurface model, boundary layer scheme, and 
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condensation (i.e. surface-atmosphere mass exchange) cycle. All of these 
features have already been included in many exoplanet models; it is merely 
a matter of integrating them into our present model. 
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Appendix: The validity of the hydrostatic approximation 

The vertical momentum equation (neglecting friction) is 



DW 2nucost- {u2 + v2) =- 1 ~^-g, (A-l) 



Dt a p dz 



where p is the atmosp 



techniques (e.g. 



leric d ensity and z is altitude. Using scale analysis 



Holtonl . 



19921 ). Eq. lA-ll may be rewritten by its scaling terms 



as 



™ Sja -**- G , (A . 2) 

L Q a Q p n 

where U is a characteristic velocity scale (taken to be 100 m s _1 ), W is a char- 
acteristic vertical velocity scale in height coordinates (taken to be 0.1 m s _1 ), 
L is a characteristic horizontal length scale (taken to be 10 7 m), f a is the av- 
erage Coriolis parameter (taken to be 10~ 4 s _1 ), a Q is the order-of-magnitude 
planetary radius (taken to be 10 7 m), P Q is the surface pressure (here we use 
the nominal case of 10 5 Pa), p Q is the representative surface density (here 
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taken to be for the 10 5 Pa case, 1 kg m 3 ), if is a representative vertical 
depth (here taken to be the scale height for the 10 5 Pa case, 10 km), and 
G is an order of magnitude for the surface gravitational acceleration (taken 
to be 10 m s~ 2 ). Note that, following the method of scale analysis, we are 
concerned with the order-of-magnitude values of the terms in Eq. \K-2\ not 
their precise values. Then, from left to right, the values of the terms in IA-21 
are: 10~ 6 , 10~ 2 , 10~ 3 , 10, and 10. This indicates that the two terms on 
the right hand side of I A- II are by far the largest terms, i.e., the hydrostatic 
balance approximation is valid. Note that the first term on the right hand 
side depends on surface pressure, which we vary in our parameter sweep. 
However, as surface pressure increases, so will p and H due to the increased 
greenhouse effect, and this term will remain the same order of magnitude. 
The hydrostatic balance approximation is then still valid. 
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Tabic 1: Selected GCM parameters 



Parameter 



Value 



Surface gravitational acceleration, g 

Orbital eccentricity 

Orbital semimajor axis 

Stellar mass 

Rotation rate (Q) 

Obliquity 

Solar constant 

Surface emissivity, e 

Radius 

Friction time scale (1/fc/) 
Specific heat at constant pressure 
Specific latent heat of vaporization 
Vaporization reference pressure (Ti) 
Vaporization reference temperature (pi) 
Specific gas constant 

Ratio of specific gas constant to specific heat 
Long wave column optical depth at 10 5 Pa (r 0l 
Radiative relaxation time scale (1//ct) 



8.93 m s- 1 


2.1426xl0 9 m 
3.1243 xlO 29 kg 
4.602xl(T 5 s- 1 
° 

21519 W m- 2 
1.0 

17059 km 
1 Earth-day 
1850 J kg- 1 K- 1 
22.6 10 5 J kg" 1 
1.01325 10 5 Pa 
373 K 
461 J kg- 1 K- 1 
1/4 
1.0 

12.6 Earth-days 
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p =10 5 Pa, A=0.1 
p s =10 5 Pa, A=0.7 
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Days after initialization 

Figure 1: Globally averaged kinetic energy per unit mass from select GCM runs. The 
kinetic energy begins from a rest state (zero) and increases as the model spins up. The level 
of kinetic energy then fluctuates due to statistical variations in the model, but maintains a 
flat trend. The results presented subsequently in this work are taken from this flat period. 
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Figure 2: Zonally averaged potential temperature (K) from GCM (global mean surface 
pressure=10 5 Pa, surface albedo=0.4 case). Contour interval is 100 K. Time averaging 
was performed over 400 Earth-days. 
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Figure 3: Temperature profiles from GCM (global mean surface pressure=f0 5 Pa, surface 
albedo=0.4 case). The profiles are time averaged over 400 Earth-days. Two latitudes are 
shown: the equator and 60°(midlatitudcs). Local times (LT) are plotted for the nightsidc 
(0 LT), dayside (12 LT), and morning (6 LT) and evening (18 LT) terminators. 
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Figure 4: Horizontal winds (arrows, in m s _1 ) and vertical winds in pressure coordinates 
(colored contours, in increments of 10~ 2 Pa s) from GCM simulation with global mean 
surface pressure=10 5 Pa and surface albedo=0.4. Time averaging was performed over 400 
Earth-days. The winds are shown as an example of the circulation in the free atmosphere 
at a level 0.068 (pressure 6800 Pa). Note that in pressure coordinates positive vertical 
wind is downward. 
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Figure 5: Zonally averaged zonal winds (m s _1 ) from GCM (global mean surface 
pressure=10 5 Pa, surface albedo=0.4 case). Contour interval is 25 m s _1 . Time aver- 
aging was performed over 400 Earth-days. The westerly equatorial jet has a width of 
about 30-45°. 
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Figure 6: Latitudinal transport of eastward eddy momentum per unit mass (m 2 s~ 2 ) from 
GCM (global mean surface pressure=10 5 Pa, surface albedo=0.4 case). Contour interval 
is 250 m 2 s~ 2 . Time averaging was performed over 400 Earth-days. 
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Figure 7: Zonally averaged mass stream function (10 12 kg s _1 ) from GCM (global mean 
surface prcssurc=10 5 Pa, surface albedo=0.4 case). Contour interval is 0.1 x 10 12 kg s _1 . 
Time averaging was performed over 400 Earth-days. Positive flow is clockwise. 
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Figure 8: Temperature profiles as a function of global mean surface pressure and surface 
albedo. GCM refers to the dynamical temperatures output from the GCM, and EQ refers 
to the radiative-convective equilibrium temperature. The profiles are time averaged over 
400 Earth-days. All profiles are shown for the equator and 12 hours LT. 
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Figure 9: GCM temperature (K) as a function of global mean surface pressure: (a) p s = 
10 4 Pa, (b) p s = 10 5 Pa, (c) Ps = 10 6 Pa, and (d) p s = 10 7 Pa. All cases have A = 0.4, 
are time averaged over 400 Earth-days, and are shown at the same a level 0.068. Contour 
interval is 5 K. 
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Figure 10: GCM temperature (K) as a function of surface albedo: (a) A = 0.1, (b) A = 0.7. 
All cases have p s — 10 5 Pa, are time averaged over 400 Earth-days, and are shown at the 
same a level 0.068. Contour interval is 5 K. Taken with panel (b) of Fig. |9l these plots 
form a progression of surface albedos from 0.1-0.7. 
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Figure 11: GCM zonal winds (m s x ) as a function of global mean surface pressure: (a) 
p s = 10 4 Pa, (b) p s = 10 5 Pa, (c) p s = 10 6 Pa, and (d) p s = 10 7 Pa. All cases have 
A = 0.4, are time averaged over 400 Earth-days, and are shown at the same a level 0.068. 
Contour interval is 50 m s _1 . Positive corresponds to westerly flow. 
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Figure 12: GCM zonal winds (m s _1 ) as a function of surface albedo: (a) A = 0.1, (b) 
A = 0.7. All cases have p s = 10 5 Pa, are time averaged over 400 Earth-days, and are shown 
at the same a level 0.068. Contour interval is 50 m s _1 . Taken with panel (b) of Fig. [11] 
these plots form a progression of surface albedos from 0.1-0.7. Positive corresponds to 
westerly flow. 
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Figure 13: GCM meridional winds (m s 1 ) as a function of global mean surface pressure: 
(a) p s = 10 4 Pa, (b) Ps = 10 5 Pa, (c) p s = 10 6 Pa, and (d) p s = 10 7 Pa. All cases have 
A = 0.4, are time averaged over 400 Earth-days, and are shown at the same a level 0.068. 
Contour interval is 25 m s _1 . Positive corresponds to southerly flow. 
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Figure 14: GCM meridional winds (m s _1 ) as a function of surface albedo: (a) A = 0.1, 
(b) A = 0.7. All cases have p s = 10 5 Pa, are time averaged over 400 Earth-days, and 
are shown at the same a level 0.068. Contour interval is 25 m s" 1 . Taken with panel 
(b) of Fig. [T31 these plots form a progression of surface albedos from 0.1-0.7. Positive 
corresponds to southerly flow. 
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Figure 15: GCM surface pressure (Pa) variation as a function of global mean surface 
pressure: (a) p s = 10 4 Pa, (b) Ps = 10 5 Pa, (c) p s = 10 6 Pa, and (d) p s = 10 7 Pa. Each 
panel is normalized by the global mean surface pressure for that case. Contour interval 
is 1% of the global mean surface pressure. All cases have A = 0.4 and are time averaged 
over 400 Earth-days. 
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Figure 16: GCM surface pressure variation as a function of surface albedo: (a) A = 0.1, 
(b) A — 0.7. All cases have p s = 10 5 Pa and are time averaged over 400 Earth-days. Each 
panel is normalized by the global mean surface pressure for that case. Contour interval is 
1% of the global mean surface pressure. Taken with panel (b) of Fig. [T5l these plots form 
a progression of surface albedos from 0.1-0.7. 
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Figure 17: Model H2O transmission spectra calculated from GCM pressure-temperature 
output. The units for p s are in bars (1 bar = 10 5 Pa). 
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